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Abstract 

The p + 6 Li — ► n + 7 Be reaction has been investigated with an emphasis on the n meson 
and 7 Be interaction in the final state. Considering the 6 Li and 7 Be nuclei to be a-d and a- 3 He 
clusters respectively, the reaction is modelled to proceed via the p + d[a] — > 3 He[a] + n reaction 
with the a remaining a spectator. The n meson interacts with 7 Be via multiple scatterings on the 
3 He and a clusters inside 7 Be. The individual r/- 3 He and n-a scatterings are evaluated using few 
body equations for the 77-3N and 7/-4N systems with a coupled channel 77-N interaction as an input. 
Calculations including four low-lying states of 7 Be lead to a double hump structure in the total 
cross section corresponding to the L = 1 ( J = (1/2)^,(3/2)") and L = 3 ( J = (5/2)^,(7/2)") 
angular momentum states. The humps arise due to the off-shell rescattering of the n meson on the 
7 Be nucleus in the final state. 



PACS numbers: 25.40.Ve, 21.85. +d, 25.10,+s 



I. ETA MESON INTERACTIONS WITH LIGHT NUCLEI 



The past few years have seen extensive investigations of the 77 meson producing reactions 
at close to threshold energies. The experiments are aimed either at directly searching [l| for 
the possible existence of eta-mesic nuclei as a result of the strong attractive nature of the r\- N 
interaction 2J or studying the final state eta-nuclear interaction to eventually conclude on the 
existence of eta-mesic nuclear states 3, 4J]. A common feature of the data on 77 production 
near threshold is the strong deviation of the cross sections from phase space. It can be 
understood as a manifestation of the strong attractive 77-N interaction (arising basically due 
to the proximity of the 77-N threshold to the Sn resonance N*(1535)). Experiments on 77 
production have been performed in nucleon-nucleon collisions and have been extended to 
proton-deuteron and deuteron-deuteron collisions too 5|. Though the focus of reactions 
such as the p + d— > p + d + 77 and p + d — > 3 He + r\ is on investigating possibilities 
of bound states of eta mesons and 2-3 nucleon systems, theoretical studies of the reaction 
mechanism revealed interesting features too In the p-d collisions, the production near 
threshold is found to be dominated by a two-step mechanism where the large momentum 
transfer in producing the 77 meson is shared by three nucleons. These findings naturally led 
to the curiosity of what happens when the 77 interacts with more than three or four nucleons 
in the final nucleus. With this motivation, measurements of the p + 6 Li — » rj + 7 Be 
reaction were carried out by the Turin group in 1993 [7J] at an incident energy of 683 MeV. 



A theoretical study of this reaction a 
a + 6 Li — > b + 7 Li was performed in 



ong with others of the type, a + 6 Li — > b + 7 Be and 



]. Part of the emphasis of this work was on obtaining 
the right form factors for 7 Li (and 7 Be) and the interactions of the mesons with the nuclei 
in the final states were neglected. The interest in the p + 6 Li — > rj + 7 Be reaction revived 
once again by the recent proposal of studying this reaction at COSY, Julich, at an incident 
energy of 673 MeV [9j. 

Having performed detailed theoretical studies of the p + d^p + d + r/ and p + d — > 
3 He + 77 reactions 6] and the 77 meson interactions with the deuteron, 3 He and 4 He nuclei 0, 
we now develop a model to study the interaction of rj mesons with a 3 He- 4 He cluster, namely, 
the 7 Be nucleus. The 7?- 7 Be interaction is then incorporated in a theoretical calculation of 
cross sections for the p + 6 Li — > rj + 7 Be reaction, with four possible low-lying states of the 
7 Be nucleus. An analysis of its effects on the p + 6 Li — > 77 + 7 Be cross sections presented 
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here, should be useful in motivating further experimental studies of this reaction. 



II. CLUSTER MODEL APPROACH 

Based on literature which supports considering the 7 Be nucleus as a bound state of an 
alpha ( 4 He) and 3 He [loj], we model the 7/- 7 Be final state interaction in the form of a three 
body problem of the ?7- 3 He- 4 He interaction. Regarding the 6 Li too as a cluster of an alpha 
and a deuteron the p + 6 Li — > r\ + 7 Be reaction is considered to proceed through 
the p + d — > 3 He + rj reaction with the a remaining a spectator. Besides, the present 
work focuses on the low energy region of 77 production where (a) the p + d — > 3 He + 77 
production amplitude is large (i.e. the p and d interact strongly to produce an 77) and (b) 
the cluster picture of low lying levels of 7 Be and 6 Li is reasonably good. There exists in 
principle, the possibility of considering the deuteron as a spectator. However, a reaction of 
the type p + a — > 5 Li + 77 followed by the combination of 5 Li and a deuteron to form the 
7 Be nucleus, does not agree with the cluster model approach here (since 7 Be is hardly known 
to be a cluster of d + 5 Li, where 5 Li is in fact unstable). The possibility of an intermediate 
p + a-^p + a+ rj reaction, with further steps of the final state p from this reaction 
combining with the spectator deuteron, i.e., p + d forming a 3 He which eventually combines 
with the final state a (from the above reaction) to form 7 Be still remains. However, this is 
not a practical option with no information available on the p + ct— > p + a + 77 reaction. 



A. Elastic scattering of eta mesons on 7 Be 

As mentioned above, we consider the 7 Be nucleus as a two body system made up of a 
3 He and 4 He nucleus and construct an elastic transition matrix for the three body problem 
of an 77 meson, 3 He nucleus and an a ( 4 He). The individual scattering of the r\ meson on 
3 He and a is evaluated using the t-matrices constructed earlier by the present authors jjj. 
These t-matrices are numerically evaluated using few body equations and include the off- 
shell rescattering of the 77 on the nucleons inside 3 He and 4 He. The r?- 3 He t-matrix is tested 
to reproduce the p + d — > 3 He + 77 cross section reasonably well 6j. 

To formulate the three body problem of the 77- 3 He- 4 He interaction, let us define fi and 
r*2 to be the coordinates of the 3 He and 4 He nuclei respectively, with respect to the mass- 
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7 centre of mass system. Defining the internal Jacobi coordinate of the relative distance 
between the 3 He and 4 He nuclei as x, one can see that r\ = a\x and f 2 = a 2 x with a\ = 4/7 
and a 2 = —3/7. The 7/- 7 Be t-matrix is then written as, 

T v - 7Be {k',k,z) = J d 3 x \^ 7 L (x)\ 2 [T^k', k, a x f, z) + T 2 (k',k,a 2 x,z)} (1) 

where Ti(k', k, a\x, z) and T 2 (k', k, a 2 x, z) are the medium modified t-matrices for the off- 
shell t] scattering on the bound 3 He and 4 He respectively. ^\{x) represents the cluster wave 
function of 7 Be with angular momentum L and z = E — \eq\ + ie, where, E is the total 
^-nucleus energy in the centre of mass and Eq is the energy required for the break up of 7 Be 
— > 3 He+ 4 He. The in-medium ^- 3 He and ^- 4 He t-matrices are written using a Faddeev type 
decomposition 12], namely, 



Ti(k',k,aiX,z) = ti(k',k,(nS,z) + J ^ - k^l^ ^ T 3$">^> a & z )> ( 2 ) 

where 2 = 1,2 and the indices 1 and 2 correspond to the 3 He and 4 He t-matrices respectively. 
The ti's represent the single scattering terms and are the matrices for purely elastic ?^- 3 He 
and 7]- 4 He scattering. They have the form, 

U(k', k;fl;z) = U(k', k;z)exp[i(k -k') • n], (3) 

with fj = ciiX and i = 1,2 as mentioned above. At the low energies considered here, the 
f]-N interaction is dominated by the Sn resonance N*(1535) and hence we perform a partial 
wave expansion and retain only s-waves, which reduces (J2]) to the following two equations 
for r/- 3 He (Ti) and 77- 4 He (T 2 ): 

T 1 (k',k,a 1 x, z) = tx(kf,k,aix,z) + J k" 2 t ~~j^W^' T ^ k "' k > Q2X ' z ) ( 4 ) 

T 2 (k',k,a 2 x,z) = t 2 (k' ,k,a 2 x, z) + J k" 2 T i(k", k, aix, z) , (5) 

where ti(k', k ; ; z) in Eq. ([3]) has reduced to U{k' , k, a^x, z) written in terms of the spherical 
Bessel functions j (aixk) and j (aixk'). Replacing the equation for T 2 in T 1; we obtain a 
recursive relation for 7\ as follows, 

Ti(k', k, aix, z) = t^k'.k.axx.z) + j k" 2 ^^^^^ '> k ' a<2X ' ^ ( 6 ) 

r r dk"dk //2 ~ 2 ti(k', k", aix, z) t 2 (k", k, a 2 x, z) Ti(k, k, a\x, z) 
+ J J 4vr 4 ( z - k" 2 /2jj) (z - ~k 2 /2fi) ' 



Once T\ is evaluated from (Q, it can be replaced into the equation for T 2 and the two can 
be substituted in ([!]) to evaluate the ?7- 7 Be t-matrix. Thus Q is evaluated retaining only 
s- waves in T\ and T 2 . 7\ is evaluated numerically using the 7/- 3 He and ^- 4 He t-matrices, 
t\ and £2 respectively, as inputs. £1 and t 2 are themselves evaluated numerically using few 
body equations and an input coupled channel rj-N t-matrix. Details of this formalism can 
be found in [31]. The two models of the elementary r/N t-matrix used here will be discussed 
in the next subsection. 



1. Models of elementary i]-nucleon scattering 

The coupled channel rj-N t-matrix which is required for the evaluation of the i]- 3 B.e and 
r/- 4 He t-matrices, t\ and t 2 , is taken from two different models available in literature. In 



13J, a transition matrix including the 7rN and r/N channels with the N*(1535) resonance 
playing a dominant role was constructed. It consisted of the meson - N* vertices and the 
N* propagator as given below: 

t vN ^ vN (k', k-z) = V ( z ) {k \ N f n ( 7 ) 

with, r N „ (z) = (z — M — Hn(z) — S^(-z) + ie where H a ( z ) ( a — ^1 v) are the self energy 



contributions from the irN and r]N loops. In 13j elastic and inelastic ?7-deuteron scattering 



as well as rj photoproduction on the deuteron was studied using this rjN model. We shall use 
two parameter sets available, one which leads to a scattering length of a^N = (0.75, 0.27) fm 
and another which leads to a V N = (0.88, 0.41) fm. We shall refer to this model henceforth 
as Model A. 

Model B used in the present work is taken from [l4j. The t-matrix for rjN —>■ rjN is 
written in a separable form as, 

t vN ^ vN (k', k;z) = v(k')t m (E)v(k) , (8) 

where the on-shell part, t m (E), is described in an effective range approximation as, 

t m (E)-i + iqv(q) 2 = l - + r -± q * + S q\ (9) 

a Zi 

The off-shell form factors have the form, v(k) = 1/(1 + A; 2 A 2 ), where A is the length 
parameter in this model. The parameter sets in this model are obtained from a fit to the 



ttN — > 7riV, 7riV — > 77A, 7iV — > 7rA and 7 A — > ?yiV data. The parameters required in Eqs. 
(JHJ) and (jUJ) above can be found in [141 ] . We choose four parameter sets with r/iV scattering 
lengths of (0.88,0.25) fm, (0.77, 0.25) fm, (0.51,0.26) fm and (0.4, 0.3) fm. 



2. Kinematics of the many body problem 

The choice of kinematic variables in a multiple scattering formalism where the many 
body scattering matrix is written in terms of a two body matrix is not unique in literature. 
The two most commonly used approaches are (i) the fixed scatterer approximation (FSA) 



and (ii) the on-energy shell impulse approximation (OEI) 15J (related to yet another ap 



proach, namely, the fixed impulse approximation (FIA) [16|). The difference between these 
approaches (which becomes important at intermediate and high energies) lies in the fact 
that in the FSA, the struck nucleon is assumed to recoil with the target as a whole, while 
in the OEI, it recoils freely. This means that in the FSA, the two-body operator does not 
follow from the equation for a free two-body t-matrix, but rather contains the mass of the 
nucleus (in which the two-body system is embedded) in the kinetic energy. Though this 
seems to be mathematically sound, such a two body t-matrix can have physically undesired 



features. For example, in 17] it was found from a phase shift study that such a two-body 
t-matrix may not display the resonant behaviour which it would be expected to. In the same 
work, in connection with 7r-nucleus scattering, it was found that considering a free two-body 
matrix simulated the contribution of continuum states (otherwise neglected in that work) 
and brought theory in closer agreement with data. 

Though the differences arising from the particular approaches used above may not be 
crucial at the low energies considered in the present work, the many-body problem con- 
sidered here is a bit more complicated and hence a small discussion is in order. In the 
works mentioned above, one studies the differences of the approaches involved, in a multiple 
scattering of a hadron on the individual nucleons in the nuclear target. Here however, the 
problem appears to be that of one multiple scattering problem embedded inside another. 
The 1] mesons scatter off the 3 He and 4 He nuclei as in a three body multiple scattering 
problem (of the ?7- 3 He- 4 He system). However, the individual r/- 3 He and r/- 4 He scatterings 
are represented by t-matrices for the multiple scatterings of the rj on the three and four 
nucleons inside 3 He and 4 He. We choose then to work in a framework where we start by 
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evaluating the momenta k and k! in the 7/- 7 Be centre of mass frame and then evaluate the 
individual r]- 3 B.e and r]- A B.e t-matrices at a Lorentz boosted energy-momentum in the 7/- 3 He 
and r]- 4 He centre of mass systems. The energy in the in-medium propagators (inside 7 Be) 
is however taken to be in the i]- 7 Be centre of mass system. 



3. Cluster wave functions of 6 Li and 7 Be 

Since the energy spacing between the first four low-lying levels of the 7 Be nucleus is small, 
we include the contribution from these four levels. We consider the angular momentum 
states with L = 1 and L = 3 corresponding to the J = (3/2~ , l/2~) and J = ( 7/2~ , 5/2~) 
levels respectively. The cluster wave functions for 7 Be and those for 6 Li required in the 
production amplitude of the p + 6 Li — > rj + 7 Be reaction are (1) generated using a Wood- 
Saxon potential [18( and (2) taken from a Green's function Monte Carlo (GFMC) variational 
calculation with the Urbana potential [19]. In the GFMC case, the wave function for 7 Li is 
used to represent the 7 Be one. The a and deuteron cluster in 6 Li is assumed to be in the 
L = state. 



B. Production mechanism of the p + 6 Li — > rj + 7 Be reaction 

Assuming that the beam proton interacts with a loosely bound deuteron in the 6 Li nucleus 
to produce an r} meson and 3 He (with the a remaining a spectator) the production amplitude 
for the p + 6 Li — » i] + 7 Be reaction can be written in terms of that for the p + d — > 3 He + 77 
process. The off shell r] meson thus produced then rescatters on 7 Be (i.e. the 3 He- 4 He 
cluster). This production of the rj and its final state interaction (FSI) with 7 Be is represented 
schematically in Fig. 1, and the corresponding transition matrix is written as, 

T = (k v ,m 7 \ T p 6 Li _+ n 7 Be \ k p ] m p , m 6 ) (10) 
r dq (k v , m 7 \T v 7 Be ^ v 7 Be \q, m' 7 ) -> 
+ If J (2^F E(k v ) — E{q) + ie ( 9 ' ™ 7 1 ' ' m *> ' m ^ 

where k p and k v are the initial and final momenta in the centre of mass system. m p , nig and 
rrij are the spin projections of the proton, 6 Li and 7 Be respectively. The production matrix 
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-J7 — ~^ 'W 7^1 — V" 17 \ ^ 7*7 N H 71 ^~ 

i t ; i I 7 r p i iTii i 'ggi i 'rp) 



d ■ 'He _ _ d f 3 He 



(a) (b) 

FIG. 1: Cluster model for 77 production in the p + 6 Li — > 77 + 7 Be reaction. Diagram (a) 
corresponds to the direct on-shell 77 production and (b) to an 77 which is first produced off-shell 
and rescatters via the r/ + 7 Be — > 77 + Be process to become on-shell. 

for a relative angular momentum L between the 3 He and a is written as, 

(k v (^, m' 7 \T pGLi ^ ri 7 Be \k p ; m p , m 6 ) = j^^-j= J P x 2 dP 1 dn Pl ^ 6 (P 1 ) (11) 

£ ( J, m' 7 1 i fi, L, M) vryf(P 2 ) F L * M (A) 

—3 -> ^ 1 1 ^ 1 ^ ^ 

x ( k v (q), —k v (q)+P 2 , -, fi\ T pd ^ v 3 He | -, m p , 1, m 6 , --fc p + Pi ) , 

where, k p is the momentum of the incoming proton and —[2/6]k p + P\ of the deuteron in 
6 Li. fi is the spin projection of 3 He. The on-shell 77 meson momentum is denoted as k v 
and the off-shell one as q. [—3/7]k v + P2 or [—3/7)5 + P2 is hence the momentum of on- or 
off-shell 3 He in 7 Be. Since the a particle remains a spectator, its momentum in 6 Li and 7 Be 
is required to be the same and — [4/6] /c p — Pi = — [A/7]k v (q) — P2 (for on-shell (off-shell) rj 
production). Thus P2 = Pi + [2/3] k p — [A/7]k v (q) (where Pi and P2 are the Fermi momenta 
inside 6 Li and 7 Be respectively). The integration in ( fTTi) should in principle have been over 
both the Fermi momenta, Pi and P2, however, the above relation between them renders the 
integration over P 2 in ( TTTj) redundant. Further, when one evaluates the unpolarized cross 
sections, one sums over the spins in the final state and averages over those in the initial 
state. As a result, in such a calculation, some sums in (TXOT) and ( JTT1) become redundant. 

The T— matrix for the process, ( | T p d^n 3 He I ), is written in a two-step model from our 
earlier work fj. Considering the complexity of the present calculations which include the 
off-shell rescattering as given by the second term in (jTUI) with the 7/- 7 Be FSI and the fact 
that the two-step model of the p + d — ► 3 He + 77 is itself quite involved, we neglect the 
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effect of Fermi motion on ( | Tpd^n^He I ) and hence take it out of the integral over P\ in (fTTI) . 
The momentum space wave functions are expressed in terms of Fourier transforms of their 
radial forms. Thus the integral in momentum space is transformed to that in coordinate 
space. All this simplifies Eq. ffTTl) to a good extent and it can be written as, 



T p e Li ^ v r Be \...) = i^V^Y, Y LM{Q)FL{Q){J,m' 7 \\, fl, L, M) 
X ( ••• I Tpd-^^He | ••• ) 



where, 



F L (Q) = / r 2 dr^l 7 (r)j L (Qr)^ (r) 
Jo 

is the transition form factor for 6 Li — > 7 Be, with momentum transfer Q 



(12) 
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FIG. 2: The 6 Li- Be transition form factor for angular momentum states with L = 1 and L = 3. 
The solid line corresponds to the GFMC variational wave functions with the Urbana potential 
for L = 1 and the dashed, dot-dashed and double dot-dashed lines are those generated using a 
Woods-Saxon potential (for L = 1 and L = 3). 

for example in the off-shell case. Though not written explicitly, the transition form factor 
depends on the total angular momentum J since the radial wave function in 7 Be depends, 
even if mildly, on J. In Fig. 2 we present the two form factors with L = 1 (J = l/2~, 3/2~) 
and L = 3 ( J = 5/2~, 7/2~) required in the present work, using two different prescriptions 
of the nuclear wave functions mentioned in the previous section. With the L = 1 levels, 
J = l/2~ and J = 3/2~ being very close to each other, the difference between the two 



L = 1 form factors is not visible in the figure. The form factor enters (fTUI) via (TT2|) as (fTT|) is 
simplified to (fT2l due to the neglect of the Fermi motion in the p + d — > 3 He + rj t-matrix. 
It is evaluated at Q = ~k v — |/c p in the first term of (1101) whereas over a range of momenta 
in the second term of ( fTUj) where it appears inside the integral. The dotted vertical lines in 
Fig. 2 indicate the relevant range corresponding to the beam energies of the present work. 
We shall see later how this difference between the two types of form factors affects the cross 
sections. 

III. RESULTS AND DISCUSSIONS 

In what follows, we shall present the cross section calculations for the p + 6 Li — > 77 + 7 Be 
reaction with an emphasis on the i]- 7 Be FSI. The transition matrices for the elementary 
processes rr + N — * rj + N and rj + N — * 77 + N which enter as inputs to the production 
and FSI matrices are chosen from coupled channels calculations. The calculations in Figs 
3-7 are done within Model A with a choice of parameters which corresponds to a scattering 
length of &r[N — 0.88 + 2 0.41 fm. Within the two-step model of the p + d — > 3 He + 77 
reaction we use here, the data on this reaction were reproduced well with this choice of the 
scattering length. 

In Fig. 8, a comparison of the total (inclusive) cross sections within models A and B of 
the rjN — > i]N t-matrix and using different sets of scattering lengths is made. 

A. The p + 6 Li -> 77 + 7 Be cross sections for different 'Be levels 

As mentioned earlier, the p + 6 Li — > 77 + 7 Be reaction is studied with four possible final 
states of the 7 Be nucleus. The proton beam energies are chosen to study the cross sections 
up to excess energies of about 20 MeV above threshold. Due to the differences between the 
masses of the various 7 Be levels, the threshold for the reaction corresponding to each level 
differs. In Table I, as an example, we list some of the beam energies at which we evaluate 
the cross sections, along with the corresponding excess energies in order to facilitate the 
understanding of the plots later. 

In Figs 3 and 4, the angle integrated total cross sections for the p + 6 Li — > 77 + 7 Be 
reaction are shown as a function of the proton beam energy. The dashed lines are the cross 
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FIG. 3: Total cross sections as a function of the proton beam energy for different states of the 7 Be 
nucleus. The dashed lines are plots without the inclusion of the ?7- 7 Be final state interaction. The 
wave functions for 6 Li and 7 Be are generated using the Woods-Saxon potential. 
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FIG. 4: Same as Fig. 3 except for the fact that the wave functions for 6 Li and 7 Be are from the 
GFMC variational method with the Urbana potential. 

sections evaluated using only Fig. 1(a) corresponding to the first term in Eq. ( jTUl) . As is 
evident from these plots, inclusion of Fig. 1(b), i.e., the rescattering of the rj meson with 
the 7 Be nucleus (or in other words the final state rj- 7 Be interaction) drastically affects the 
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TABLE I: Beam and excess energies for different levels of the Be nucleus 



Beam energy 
(MeV) 


Excess energy = E^_7^ e - M^Be - (MeV) 


L = 1 


L = 3 


J = (3/2)- 
ground state 


J = (1/2)" 


J = (7/2)" 


J = (5/2)- 


658.8 


0.0003 








659.8 


0.791 


0.361 






663.8 


3.954 


3.524 






669 


8.056 


7.626 


3.426 


1.376 


673.8 


11.855 


11.425 


7.225 


5.175 


683 


19.11 


18.679 


14.478 


12.428 



shape and magnitude of the cross sections near threshold. 
1. The off-shell rescattering contribution 

The second term in Eq. (TTTJT) which corresponds to rescattering (Fig. 1(b)), consists 
in principle of the scattering of on- as well as off-shell i] mesons on 7 Be. As is generally 
expected at low energies, Fig. 5 clearly shows that neither the plane wave scattering (Fig. 
la) nor the pole term in Eq. (|1~0|) (the on-shell 7/- 7 Be rescattering) is responsible for the near 
threshold cross section hump. It is the principal value of the integral in the second term which 
corresponds to the scattering of off-shell eta mesons produced in the p + 6 Li — > rj + 7 Be 
process on 7 Be that gives rise to this hump. Though the dominance of the principal value 
is expected, a hump-like structure in the total cross section, due to final state 7/-nucleus 
interaction is not so natural to expect. For example, such an effect was not observed in our 
previous studies of the pd — > pdr] and pd — > 3 He rj reactions. 

In order to understand this phenomenon, we re-write the principal value integral in 
the second term of Eq. ( jTTJl) (corresponding to off-shell rescattering) in a split form. Let 
us express the integral / dq = J dQ q / °° q 2 dq as / dQ q J qp,dE(q) where, q 2 = 2fiE(q). 
Dropping the spin projection dependence in the notation for convenience, we define the 
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Off-shell rescattering 

On-shell rescattering 

7 

without T|- Be interaction 



665 670 675 680 
Beam energy (MeV) 



FIG. 5: Contributions of the plane wave (dashed line), on-shell rescattering (dotted line) and off- 
shell rescattering (dot-dashed line) terms to the total cross section (solid line) for the p + 6 Li — > 
T] + 7 Be reaction with 7 Be in its ground state (J = 3/2"). The cluster model nuclear wave functions 
are generated using the Woods-Saxon potential. 

rescattering term alone (second term in (TTUl) ) as a function, G(k p ,k v ) plus the pole term, 
such that, 

dq (k v \T v 7 Be \q) 



(2tt) 3 E(k v ) - E(q) + ie 
= v[dQ q F{k p , k v , q) [ dE{q) 



( q\ T p 6 Li ^ v 7 Be | k p ) 
f(E(q)) 



(14) 



pole term . 



E{k v ) - E(q) 

where the functions, F and / summarize the angle and energy dependence (except the 
energy dependence in the denominator) respectively of the vector q in the integrand. We 
now break up part of the principal value integral into three parts, considering the region till 
and after the pole value as follows: 



E(K)-& f(E(q)) 



+ 



E(k v ) - E(q) 
f(E(q)) 



dE{q) + 



2E 



(*,) f(E(q)) 



'E(k v )+s E(k v ) - E(q) 
dE{q) = h + I 2 + J 3 . 



dE{q) 



(15) 



hE{k v ) E(k ri ) - E(q) 

In order to get some idea of the behaviour of the complicated integrals in (fT4l and (|T5l) which 
we perform numerically, we try to analyse these integrals analytically in a very simplistic 
way. First, we perform an expansion of f(E(q)), such that, 

f(E(q)) ~ f(E(k v )) + (E(k v ) - E(q))f'(E(k v )) + ... 
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and hence, 

f(E(q)) ^ f(E(k v )) 



+ f'(E(k v )) + 



E(k v ) - E(q) E(k v ) - E(q) 
We further assume that only the first term on the right in the above expansion is important. 
Not always can one expect the energy dependence of f(E(q)) to be such that we can write 
it with only a constant f(E(k v )) appearing in the first term as above. If the pole is however 
close to zero, then the energy region — > 2 E(k v ) is a small region around the pole where 
f(E(q)) may vary little and the expansion made above becomes a reasonable approximation. 
Since f(E(k r) )) is a constant, when we integrate, we get, 



E(k v )-S 



+ ]n(E(k n ) - E(q)) 



2E(k v ) 

+ h, (16) 



where obviously the first two terms cancel. The integral I in (JTSJ) will then be dominated only 
by ^3 whose behaviour will decide the shape of the off-shell rescattering term. I3 depends on 
the form of the function f(E(q)), which is a product of the p + 6 Li — > 77 + 7 Be production 
t-matrix and the r/- 7 Be elastic t-matrix. The latter is peaked at small k v due to the fact 
that the 7 Be wave function contributes at large r. The -Z3 part of the off-shell rescattering 
term which grows as one approaches energies close to threshold gives rise to the hump due 
to FSI in the total cross section. The above analysis will not be valid when one is far away 
from threshold. 



2. Single scattering versus multiple scattering to all orders 

In Fig. 6, for the sake of completeness is shown the cross section evaluated using only the 
single scattering terms in the ?7- 7 Be t-matrix (i.e. the first terms in Eqs (jU-EJ)) as compared 
to that using the full coupled channel t-matrix which includes the multiple scattering to all 
orders (Eqs (JH-E])). The dot-dashed line is thus the calculation corresponding to an r/- 7 Be 
interaction where the rj meson scatters once on each of the 3 He and 4 He nuclei inside 7 Be 
and then proceeds to become on-shell in the final state of the p + 6 Li —>■ rj + 7 Be reaction. 
The solid line is due to the ?7- 7 Be interaction which involves multiple scattering of the 77 on 
the two clusters, 3 He and 4 He, and its propagation inside 7 Be in between these scatterings. 
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FIG. 6: Total cross sections evaluated using only the single scattering term (dot-dashed) in the 
r/- 7 Be interaction, the full calculation (solid line) involving 77 meson rescattering and propagation 
inside 7 Be to all orders and the one without the inclusion of the ?7- 7 Be interaction (dashed line). 
The calculation is done for the ground state of 7 Be (J = 3/2~). 

3. Angular distributions 

In Fig. 7, we plot the angular distributions for this reaction for the ground state (J = 
3/2~) of 7 Be. The angular distributions are nearly isotropic as has always been found in the 
rj producing reactions near threshold. The FSI between the r\ meson and 7 Be is responsible 
for raising the magnitudes of the cross sections only close to threshold. At an excess energy 
of 19 MeV, the effect of the FSI is highly reduced and the theoretical prediction in Fig. 7 
is around 0.2 nb/sr as compared to the Turin data point of 4.6 ± 3.8 nb/sr. Even if we add 
up the cross sections corresponding to other values of J, i.e., J = 1/2", 5/2" and 7/2", the 
theoretical prediction with FSI is 0.51 nb/sr. Though not very clear at the moment, one 
could speculate that the deficit in the theoretical prediction as compared to data is probably 
due to some missing reaction mechanism not included in the present cluster model approach. 
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FIG. 7: Angular distributions for the p + 6 Li — > r\ + 7 Be reaction at different excess energies 
(corresponding to different beam energies). The data point is by the Turin group 7j. Solid and 
dashed lines correspond to the calculations with and without the inclusion of the ry- 7 Be interac- 
tion respectively. The calculation is done for the Be ground state and using Woods-Saxon wave 
functions. 



B. Inclusive cross section 



Having analysed and deciphered the origin of the hump in the total cross sections near 
threshold, we now plot in Fig. 8, the sum of the total cross sections, YLj v'l+^^^ + i^i 
where the superscript J refers to the state in which the 7 Be in the final state is produced. 
The four different levels of 7 Be considered here have a different value of J and also a different 
mass. In this sense, the summed cross section represents an inclusive cross section of the 
p + 6 Li — > 7] + X reaction with X being any of the low-lying states of 7 Be. We have already 
seen that within Model A, the close to threshold ?7- 7 Be interaction causes a sharp rise (the 
hump) in the cross section. Now, in the case of the L = 3 states, as can be seen from Table 
I, the threshold is shifted by about 10 MeV as compared to the L = 1 states. As a result, 
the L = 3 humps appear and contribute at a higher beam energy. Since the hump is a close 
to threshold phenomenon, the L = 1 cross sections already fall down a lot before the L = 3 
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FIG. 8: Total cross sections for the inclusive reaction p + 6 Li — > r\ + 7 Be. Inclusive implies that the 
cross section for four J values of 7 Be, namely, l/2~ , 3/2~, 5/2~ and 7/2~ have been summed. The 
upper and lower panel display results within two different models of the r]N interaction used for 
evaluating the ry- 7 Be FSI as described in the text. The dotted line without any humps corresponds 
to the calculation without the r/- 7 Be final state interaction. 

humps begin. The result is a double hump structure in the inclusive p + 6 Li — > rj + 7 Be 
reaction. 

The cross sections evaluated in Model B, have a similar form as those in Model A except 
for two values of the scattering lengths, namely, (0.77,0.25) fm and (0.4,0.3) fm. For these 
two values, the cross sections do not display a very prominent second hump. After a careful 
reading of Ref. [ijj] (which lists the parameter sets for Model B, lower panel in Fig. 8) 
one notices that the parameter sets for these two particular curves have been obtained after 
removing the jN — > rjN data from the analysis. Thus these two curves are different from 
all the rest in Fig. 8, in the sense that these parameter sets are not constrained by the 
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photoproduction data. Model A, however, was used in [13] to calculate the cross sections 
for the jd —>■ rjd and 7c? — > 7]X reactions. Hence, Model A and the two curves with double 
humps in Model B are constrained by the eta photoproduction data which is probably 
important for the behaviour of the rjN t-matrix at energies away from threshold. 



C. Comment on possible eta-mesic 7 Be states 

Finally, before ending our discussion of the results, in Table II, we give the rj- 7 Be scattering 
lengths corresponding to different values of the rjN scattering lengths in models A and B. 
The r/- 7 Be scattering length is evaluated from the elastic rj- 7 Be t-matrix at zero energy as 
follows: 

/iTV 7Be (0,0,0) 



dr]-7Be 



2tt 



(17) 



where \i is the reduced mass of the r]- 7 Be system. The smallest scattering length a V N used 



TABLE II: ry- 7 Be scattering lengths, a^-jEe, for the ground state of 7 Be, corresponding to different 
values of a^jv- 



Model A 


a vN (fm) 


ary-7.Be (fm) 


0.75 + i0.27 
0.88 + z0.41 


-10.09 + *8.19 
-9.18 + i8.53 


Model B 


0.88 + i0.25 
0.77 + i0.25 
0.51 + i0.26 
0.4 + i0.3 


-20.43 + ?5.43 
- 14.52 + »14.77 
-2.03 + ill. 29 
0.29 + i6.43 



here, leads to a positive real part of the ?7- 7 Be scattering length, a v ^ 7 B e - However, this choice 
of parameters (model F in 14() is mentioned as an unconventional solution obtained after 
dropping the photoproduction data from the fits. 

The first four entries in Table II corresponding to large a V N, display a J? _ 7 B e to be large 
with negative real parts. Apart from the commonly known condition that the real part of 
the scattering length should be negative [2^], in the third reference in the authors found 
that the condition for the existence of a bound state is that lad < |a#| for an eta- nucleus 
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scattering length of (or + iai). With these two conditions it seems that the first few 
entries in Table II for the large 7]N scattering lengths support the possibility of eta-mesic 



states. In 



2l| . while investigating the connection between the ?7- 3 He scattering lengths and 



the corresponding binding energies and widths, the authors also use the above conditions 
but mention that in reality none of the above can be taken as a sufficient condition. The 
bottomline is then that it would indeed be premature to base the conclusions only on the 
signs or magnitudes of the scattering lengths. One should rather perform a better analysis 
'or the search of 77- 7 Be states using the present r?- 7 Be model and a time delay analysis as in 



3, 



22j or a K-matrix analysis as in 23j before drawing definite conclusions. 



IV. SUMMARY 



The present work aimed at performing a thorough investigation of the effects of the 
r]- 7 Be interaction in the p + 6 Li — ^ 77 -h 7 Be reaction near threshold. The work was 
partly motivated by the recent revival of interest in this reaction by the COSY-GEM 
collaboration [9J after the first measurement in 1993. This work also comes as a sequel to 
our various earlier studies on 7] meson production in light nuclei. A two-step model for the 
p + d — > 3 He + 7] reaction including the 77- 3 He interaction was tested earlier with data by 
the present authors p]. This model is used as an input to develop a cluster model approach 
for the p + 6 Li — > 77 + 7 Be reaction near threshold. The ??- 7 Be interaction is included in a 
multiple scattering formalism for an 77 scattering on a 3 He- 4 He cluster inside 7 Be. The ??- 3 He 
and - 4 He scatterings are themselves included using few body equations. The calculations 
are done for four low-lying levels of 7 Be. To the best of our knowledge, this is the most 
detailed study of the 77- 7 Be interaction in the p + 6 Li —>■ rj + 7 Be reaction performed so 
far. The interesting two hump structure in the summed total cross section (summed over 
J) hints toward a very strong near threshold effect of the 77- 7 Be interaction (especially 
of the off-shell rescattering of the rj on 7 Be) which is worth verifying experimentally in future. 

ACKNOWLEDGEMENTS 

The authors gratefully acknowledge V. Jha for the help in computing the cluster wave 
functions using the Woods-Saxon potential and K. P. Khemchandani for her help with the 
pd — > r/ 3 He programs and many useful discussions. The authors are also thankful to the 



19 



anonymous referee for his/her constructive criticism. The work done by two of the authors, 
NJU and BKJ, was supported by the Ramanna fellowship, awarded by the Department of 
Science and Technology, Government of India, to BKJ. 



[1] M. Pfeiffer et al, Phys. Rev. Lett. 92, 252001 (2004); G. A. Sokol et al, Part. Nucl. Lett. 
102, 71 (2000). 

[2] R. S. Bhalerao and L. C. Liu, Phys. Rev. Lett. 54 (1985) 865; L.C. Liu and Q. Haider, Phys. 

Rev. C 34, 1845 (1986); Q. Haider and L.C. Liu, Phys. Rev. C 66, 045208 (2002). 
[3] N. G. Kelkar, K. P. Khemchandani and B. K. Jain, J. Phys. GrNucl. Part. Phys. 32, L19 

(2006); N. G. Kelkar, K. P. Khemchandani and B. K. Jain, J. Phys. G 32, 1157 (2006); N. 

G. Kelkar, Phys. Rev. Lett. 99, 210403 (2007). 
[4] S. Wycech, Anthony M. Green and J.A. Niskanen, Phys. Rev. C 52, 544 (1995); V.A. Trya- 

suchev, Phys. Atom. Nucl. 60, 186 (1997) (Yad. Fiz. 60, 245 (1997)); C.Y. Song, X.H. Zhong 

and L. Li, PZ. Ning, Europhys. Lett. 81, 42002 (2008); D. Jido, E.E. Kolomeitsev, H. Na- 

gahiro and S. Hirenzaki, Nucl. Phys. A 811, 158 (2008); H. Nagahiro, D. Jido and S. Hirenzaki, 

larXiv:0811.45161 

[5] J. Berger et al, Phys. Rev. Lett. 61, 919 (1988); B. Mayer et al, Phys. Rev. C 53, 2068 
(1996); A. Wroriska et al, Int. J. Mod. Phys. A 20, 640 (2005); H., -H., Adam et al, Phys. 
Rev. C 75, 014004 (2007). 

[6] K.P Khemchandani, N.G. Kelkar and B.K. Jain, Nucl. Phys. A 708, 312 (2002); ibid, Phys. 
Rev. C 68, 064610 (2003); N.J. Upadhyay, K.P. Khemchandani, B.K. Jain and N.G. Kelkar, 
Phys. Rev. C 75, 054002 (2007); K. P. Khemchandani, N. G. Kelkar and B. K. Jain, Phys. 
Rev. C 76, 069801 (2007). 

[7] Scomparin E. et al, J. Phys. G: Nucl. Part. Phys. 19 L 51, (1993). 

[8] J. S. Al-Khalili, M. B. Barbara and C. Wilkin, J. Phys. G 19, 403 (1993). 

[9] M. Ulicny et al, AIP Conf. Proc. 603, 543 (2001); H. Machner, Acta Phys. Slov. 56, 227 
(2006); e-Print: |nucl-ex/0511034 
[10] T. Kajino, T. Matsuse and A. Arima, Nucl. Phys. A 413, 323 (1984). 

[11] J. V. Noble, Phys. Rev. C 9, 1209 (1974); A. C. Merchant and N. Rowley, Phys. Lett. B 150, 



20 



35 (1985). 

[12] S. A. Rakityansky et al, Phys. Rev. C 53, 2043 (1996); S. A. Rakityansky et al, Few Body 

Syst. Suppl. 9, 227 (1995). 
[13] A. Fix and H. Arenhovel, Eur. Phys. J A 9, 119 (2000). 
[14] A. M. Green and S. Wycech, Phys. Rev. C 71, 014001 (2005). 
[15] E. Kujawski and E. Lambert, Annals of Physics 81, 591 (1973). 
[16] R. Crespo, A. M. Moro and I. J. Thompson, Nucl. Phys. A 771, 26 (2006). 
[17] V. B. Belyaev, S. A. Rakityansky and J. Wrzecionko, Nucl. Phys. A 368, 394 (1981). 
[18] B. Buck and A. C. Merchant, J. Phys. G 14 L211 (1988); V. I. Kukulin, V. G. Neudatchin 

and Yu. F. Smirnov, Nucl. Phys. A 245 (1975) 429. 
[19] J. L. Forest et al, Phys. Rev. C 54, 646 (1996). 

[20] C. J. Joachain 1975 Quantum Collision Theory (North-Holland), p. 288. 

[21] A. Sibirtsev, J. Haidenbauer, J. A. Niskanen and Ulf-G. Meissner, Phys. Rev. C 70, 047001 



(2004); see also J. A. Niskanen, arXiv: nucl-th/0508021 



[22] N. G. Kelkar, M. Nowakowski, K. P. Khemchandani and S. R. Jain, Nucl. Phys. A730, 121 
(2004); N. G. Kelkar, M. Nowakowski and K. P. Khemchandani, Nucl. Phys. A724, 357 
(2003); ibid, Mod. Phys. Lett. A 19, 2001 (2004). 

[23] S. Wycech and A. M. Green, Int. J. Mod. Phys. A 20, 637 (2005). 



21 



